QCD quark cyclobutadiene and light tetraquark spectra * 



(N 

o 

(N 
-i— > 
O 

O 



J2' 

1). 
43 ■ 



>: 

oo. 



o 

(N- 



X: 



Chengrong Deng a , Jialun Ping 6 , Hui Wang 6 , Ping Zhou", Fan Wang c 
1 School of Mathematics and Physics, Chongqing Jiaotong University, Chongqing 400074, P.R. China 
b Department of Physics, Nanjing Normal University, Nanjing 210097, P.R. China and 
c Department of Physics, Nanjing University, Nanjing 210093, P.R. China 

The QCD quark cyclobutadiene (ring-like), a new color structure of tetraquark system, is proposed 
and studied in the flux-tube model with a multi-body confinement interaction. Numerical calcula- 
tions show that the light tetraquark systems (u, d, s only) with cyclobutadiene, diquark-antidiquark 
flux-tube structures have similar energies and they can be regarded as QCD isomeric compounds. 
The energies of some tetraquark states are close to the masses of some excited mesons and so in 
the study of these mesons, the tetraquark components should be taken into account. There are also 
some meson states, a, rc(800), /o(980), /o(1500), 7Ti(1400), 7Ti(1600), / 2 (1430) and if* (1410), where 
tetraquark components might be dominant. The meson states with exotic quantum numbers are 
studied as the tetraquark states. The multi-body confinement interaction reduces the energy of the 
tetraquark state in comparison with the usual additive two body confinement interaction model. 

PACS numbers: 14.20.Pt, 12.40.-y 



I. INTRODUCTION 

In the constituent quark model (CQM), mesons are as- 
sumed to be composed of qq. Although various properties 
of light mesons have been explained within this qq mini- 
mum Fock space, there are still properties of some meson 
states can not be described well by this quark model [J- 
3]. In fact, mesons might be more complicated objects 
with higher Fock components other than the lowest qq. 
The wavefunction of zero baryon number (B=0) hadron, 
if the gluon degree of freedom is neglected, should be in 
general as 
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where n = 1, 2, 3, ■ • The recent studies on meson spec- 
troscopy called for unquench the quark model, i.e., the qq 
and q 2 q 2 mixing 0-0] . Furthermore the introduction of 
tetraquark states q 2 q 2 is indispensable for the states with 
exotic quantum numbers |7H11|. In recent years, compre- 
hensive researches on tetraquark states have been carried 
out by many authors [12h21| . because Belle, BaBar and 
other experimental collaborations have observed many 
open and hidden charmed hadrons, which are difficult 
to be fitted into the conventional meson cc spectra [22j . 
The states with quantum numbers J PC = 0~~, even" 1 
and odd had been theoretically studied as tetraquark 
states [23H25j . Experimental evidences of the exotic 
states with quantum numbers J 



PC 



= V 



were accu- 



itn qui 

mulated |7Hll|. The investigations of multiquark states 
with flux-tube structures will provide important low en- 
ergy quantum chromodynamics (QCD) information, such 
as qqq and qqq interactions [26] , which is absent in ordi- 
nary hadrons due to their unique flux-tube structure. 



QCD is widely accepted as the fundamental theory 
of the strong interaction, in which color confinement is 
a long-distance behavior whose understanding continues 
to be a challenge for theoretical physics. Lattice QCD 
(LQCD) allows us to investigate the confinement phe- 
nomenon in a nonperturbative framework and its calcu- 
lations on mesons, baryons, tetraquark and pentaquark 
states reveal flux-tube or string like structure [27H30j . 
Such flux-tube like structures lead to a "phenomenolog- 
ical" understanding of color confinement and naturally 
to a linear confinement potential in qq and q 3 quark sys- 
tems. 

It is well known that nuclear force and molecule 
force are very similar e xcep t for the length and en- 
ergy scale difference jU l32j . For multi-body systems, 
the flux-tubes in a multi-quark system should be also 
very similar to the chemical bond in the molecular sys- 
tem. Among organic compounds, the same molecular 
constituents may have different chemical bond struc- 
ture, which are named as isomeric compounds. In 
hadronic world, the multiquark states with same quark 
contents but different flux-tube structures should be sim- 
ilarly called as QCD isomeric compounds. The past 
theoretical studies on multiquark states reveal various 
flux-tube structures 1331- 4211 : hadron molecular states 



[qq]i[qq]i, [qq}i[q 3 h, [q 3 hk 3 ]i and [q 3 ]i[q 3 ]i and hid- 
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[[g 3 ]s[(7 3 ]8]i and a QCD quark benzene [g 6 ]i, et al, here 
the subscripts represent color dimensions, which should 
be mixed and affect the corresponding hadron properties 
if they really exist. 

Based on the chemical benzene and the similarity be- 
tween color flux-tubes and chemical bonds, a new flux- 
tube structure, the quark benzene, for a six-quark sys- 
tem was proposed and its possible effect on NN scat- 
tering was discussed in our previous paper [4l|. In the 
present work, a new flux-tube structure for a tetraquark 
state, which is similar to the molecular cyclobutadiene 
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and is therefore called as a QCD quark cyclobutadiene, 
is proposed. The aims of this paper are: (i) to inves- 
tigate the properties of a QCD cyclobutadiene in the 
flux-tube model, this model involves a multibody con- 
finement potential and has been successfully applied to 
multiquark systems [H, HH; (ii) to study the spectra 
of light tetraquark states with two flux-tube structures 
(diquark-antidiquark and QCD quark cyclobutadiene), 
which helps us to understand the meson states beyond a 
qq configuration and will provide a new sample to study 
the mixing of qq and qqqq. The paper is organized as fol- 
lows: four possible flux-tube structures of a tetraquark 
system are discussed in Sec. II. Section III is devoted to 
the description of the flux-tube model and the multibody 
confinement potentials of a diquark-antidiquark and a 
QCD quark cyclobutadiene structures. A brief introduc- 
tion of the construction of the wave functions and quan- 
tum numbers of a tetraquark state are given in Sec. IV. 
The numerical results and discussions are presented in 
Section V. A brief summary is given in the last section. 




FIG.1 c FIG.1 d 



FIG. 1: Four possible flux-tube structures. 



II. FLUX-TUBE STRUCTURES OF A 
TETRAQUARK STATE 

In the flux-tube picture it is assumed that the color- 
electric flux is confined to narrow, flux-tube like tubes 
joining quarks. A flux-tube starts from every quark and 
ends at an antiquark or a Y-shaped junction, where three 
flux-tubes annihilate or be created [44|. In general, a 
state with N + 1-particles can be generated by replac- 
ing a quark or an antiquark in an iV-particles state by 
a Y-shaped junction and two antiquarks or two quarks. 
According to this point of view, there are four possible 
flux-tube structures for a tetraquark system as shown 
in Fig.l, where represents the position of a quark 
(antiquark q~i) which is denoted by a solid (hollow) dot, 
Yi represents a junction where three flux-tubes meet. A 
thin line connecting a quark and a junction represents 
a fundamental flux tube, i.e. color triplet. A thick line 
connecting two junctions is for a color sextet, octet or 
others, namely a compound flux tube. The numbers on 
the flux-tubes represent the color dimensions of the corre- 
sponding flux tube. The different types of flux-tube may 
have different stiffness 45j , the detail will be discussed in 



the next section. Both the overall color singlet nature of 
a multi-quark system and the ST/(3) color coupling rule 
at each junction must be satisfied. 

The flux-tube structure (a) in Fig. I is a meson-meson 
molecule state, many newly observed exotic hadrons are 
discussed in this picture p33l - l3a |. The tetraquark states 
with the flux-tubes tructure (b) generally have high en- 
ergies due to a repulsive interaction between a quark and 
an antiquark in a color octet meson. Thus this flux- 
tube structure is often neglected in the study of multi- 
quark states. However sometimes the attraction between 
two color octet mesons will lower the energies of the sys- 
tem considerably. In the case of the flux-tube structure 



(c), called as a diquark-antidiquark structure, it has two 
possible color coupling schemes, namely [[<?<?] 3 [99] 3] i and 
[[<?<?] 6 [<?<?] g]i j the latter is expected to be a highly excited 
state and therefore the 6 diquark is usually named as a 
"bad" diquark, since the interaction between two sym- 
metric quarks (antiquarks) is repulsive, thus many au- 
thors are in favor of the 3 ( "good" ) diquark picture (36l — 

m. 

The first three flux-tube structures can be explained 
as the basic structures for a tetraquark system. The last 
structure can be generated by means of exciting two Y- 
shape junctions and three compound flux-tubes from vac- 
uum based on the second or third structures. In the con- 
stituent quark model, a quark is massive. One can sup- 
pose that the recombination of flux-tubes is faster than 
the motion of the quarks. Subsequently, the ends of four 
compound flux-tubes meet each other in turn to form 
a closed flux-tube structure, a ring-yiy3y2y4, which was 
interpreted as a pure gluon state by Isgur and Paton [44j , 
as a glueball state, and discussed in the framework of 
the dual Ginzburg-Landau theory [46j]. With quarks or 
antiquarks connecting to vertexes of the square by a fun- 
damental flux tube, this picture could be explained as 
a qqqq-ghiebal\ hybrid. According to overall color sin- 
glet and SU (3) color coupling rule, the corresponding 
compound color flux-tube dimensions (^13,^32,(^24,^41) 
have six different sets: (3,8,3,8), (6,8,6,8), (3,3,3,3), 
(8,3,8,3), (8,6,8,6) and (3,6,3,6). The flux-tubes lo- 
cate in opposite sides of the ring-yiy3y2y4 have the same 
dimension, which is similar to the symmetry with the dis- 
tribution of double bonds and single bonds in a cyclobu- 
tadiene in chemistry. We thus name the flux-tube struc- 
ture (d) as a QCD quark cyclobutadiene. Of course, the 
existence of another QCD quark cyclobutadiene is also 
allowed in which two quarks or antiquarks seat neigh- 
boring positions in the flux-tubering. Certainly, more 
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complicated configuration are permitted, including more 
Y-shaped junctions and more complex topological struc- 
tures. 



III. THE FLUX-TUBE MODEL AND 
MULTI-BODY CONFINEMENT POTENTIALS 

Recently, LQCD and nonperturbative QCD method 
have made impressive progresses on hadron properties, 
even on hadron- hadron interactions [47l - l5lj ]. However, 
QCD-inspircd CQM is still an useful tool in obtaining 
physical insight for these complicated strong interaction 
systems. CQM can offer the most complete description 
of hadron properties and is probably the most success- 
ful phenomenological model of hadron structure In 
traditional CQM, a two-body interaction proportional to 
the color charges Xi-Xj and r™-, where n — 1 or 2 and rij is 
the distance between two quarks, was introduced to phc- 
nomenologically describe quark confinement interaction. 
The traditional model can well describe the properties 
of ordinary hadrons (qqq and qq) because the flux-tube 
structures for an ordinary hadron are unique and triv- 
ial. However, the traditional model is known to be flawed 
phenomenologically because it leads to power law van der 
Waals forces between color-singlet hadrons [521154} . It is 
also flawed theoretically in that it is very implausible that 
the long-range static multi-body potential is just a sum 
of the two-body ones [55|. Many papers were devoted 
to eliminate the physically nonexisting long-distance van 
der Waals force arising from the traditional models based 
on the sum of two-body Casimir scaled potentials [56l - l59l | . 

LQCD studies show that the confinement potential of 
a multiquark state is a multibody interaction which is 
proportional to the minimum of the total length of flux- 
tubes which connects the quarks to form a multiquark 
state f27l - l3"ol ] . The naive flux-tube model is developed 
based on LQCD picture by taking into account a multi- 
body confinement potential with a harmonic interaction 
approximation, i.e., a sum of the square of the length of 
flux-tubes rather than a linear one is assumed to simplify 
the calculation [4l|, |60j. The approximation is justified 
with the following two reasons: one is that the spatial 
variations in separation of the quarks (lengths of the flux 
tube) in different hadrons do not differ significantly, so 
the difference between the two functional forms is small 
and can be absorbed in the adjustable parameter, the 
stiffness. The other is that we are using a nonrelativistic 
dynamics in the study. As was shown long ago [6l|, an 
interaction energy that varies linearly with separation be- 
tween fermions in a relativistic first order differential dy- 
namics has a wide region in which a harmonic approxima- 
tion is valid for the second order (Feynman-Gell-Mann) 
reduction of the equations of motion. Combining with 
Gaussian expansion method (GEM), the flux-tube model 
including one gluon exchange and one boson exchange in- 
teractions was successfully applied to new hadronic states 
and some interesting results were obtained [42|, |43| ■ 



Within the flux-tube picture, the flux-tubes in the 
ring-structure, see Fig. Id, are assumed to have the same 
properties as the flux-tubes in the ordinary meson or 
baryon [62| . Thus in the flux-tube model with quadratic 
confinement, the confinement potentials V c and V d for 
diquark-antidiquark and cyclobutadiene structures, have 
the following forms, respectively. 



yc 



V = k 



k [(ri - yi) 2 + (r 2 - yi ) 2 + (r 3 - y 2 ) 2 

(r 4 -y2) 2 + «d 12 (yi -y 2 ) 2 ] , (2) 

4 i 

Y ( r * - yi) 2 + Y (y« - Yj) 



i<3 



(3) 



where the means that the summation is over the ad- 
jacent junction pairs on a compound flux tube, this term 
is the energy of the flux-tube ring-yiy3y2y4. The param- 
eter k is the stiffness of an elementary flux tube, while 
kn dij is other compound flux-tube stiffness. The com- 
pound flux-tube stiffness parameter n d -. depends on the 
color dimension, dij, of the flux tube [45| . 



Co 



(4) 



where is the eigenvalue of the Casimir operator as- 
sociated with the SU(3) color representation dij on ei- 
ther end of the flux tube, namely C3 = |, Cq = ^ and 
C 8 = 3. 

For given quark positions r^, the positions of those 
junctions y,, variational parameters, can be determined 
by means of minimizing the confinement potentials V c 
and V d . To simplify the formats of V c and V d after 
obtaining the positions of the junctions y,, two set of 
canonical coordinates R^ and IZi can be introduced re- 
spectively and written as, 



(5) 




(6) 



The minimums V^ in and V^ in of the confinement poten- 
tials can be divided into three independence harmonic 
oscillators and have therefore the following forms, 



y a . = 



R 2 



- R, j 



2Kdl -K\ 



1 + 2n dl 



2{n dl + K d2 ) 
1 + 2(re ( j 1 + Kd 2 ) 



r Kdi2 

1 + 2n d2 



(7) 



(8) 
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where the parameters and used to describe 

the stiffness of two sets of opposite flux-tubes in the ring- 
yiY3y2y4 due to the symmetry, respectively. Obviously, 
the confinement potentials V£ in and V£ in are multi-body 
interactions rather than the sum of two-body interac- 
tions. 

The limit going to infinity indicates that the cor- 
responding compound flux-tube contracts to a junction 
due to the requirement of the minimum of the confine- 
ment. The limit n di . going to zero indicates the rupture 
of the corresponding compound flux-tube and then a mul- 
tiquark state decays into several color singlet hadrons. 
The flux-tube structures of a multiquark state can there- 
fore change if the Kd w is taken as an adjustable parame- 
ter. In the limit or Kd 2 going to infinity, a QCD quark 
cyclobutadiene reduces to a two-color-octet meson state 
or a diquark-antidiquark state. In the limit one of Ki and 
K2 going to infinity and the other going to zero, a QCD 
quark cyclobutadiene decays into two color mesons. In 
the limit all k^'s in Fig.l going to infinity, the last three 
flux-tube structures reduce to one structure due to all 
compound flux-tubes shrink to a junction, leaving a hub 
and spokes configuration. 

Taking into account a potential energy shift A in each 
independent harmonic oscillator, the confinement poten- 
tials V^ in and V^in have therefore the following forms 



V c = k 

v rain 



V d - k 
v min 



(R? - A) + (Rj - A) 



K d 12 



1 + Kd 12 



(Rj-A) 



1 + 2 Kdl 



+ 



2[Kdi + K d 2 ) ,j,2 

i + 2(/« dl + K d2 y V;i 



1 + 2n d2 

iRi - A) 



(9) 



(10) 



where the parameters k and A are determined by fitting 
ordinary meson spectra [43j . Carlson and Pandharipande 
also considered similar flux-tube energy shift which is 
proportional to the number of quarks N [631 ] . 

One gluon exchange and one Goldstone boson ex- 
change interactions not only are important and respon- 
sible for the mass splitting in the meson spectra but also 
are indispensable for the investigations on the multiquark 
system [42j , the details of the parts of model Hamiltonian 
can be found in our previous paper [43| . 



IV. WAVE FUNCTIONS AND DEFINITION OF 
QUANTUM NUMBERS 

The flux-tube structure specifies how the colors of 
quarks and anti-quarks are coupled to form an overall 
color singlet. It is however difficult to construct the color 
wave function of the QCD quark cyclobutadiene only us- 
ing quark degrees of freedom in the framework of the 
quark models. In order to comprehensively study a QCD 
quark cyclobutadiene, one gluon exchange and one bo- 
son exchange interactions have to be included. The color 





FIG. 2: Jacobi coordinates for a q q system 



wave function of a QCD quark cyclobutadiene is there- 
fore indispensable and approximately assumed to be the 
same as that of a diquark-antidiquark structure. In the 
framework of a diquark-antidiquark structure, three rel- 
ative motions are shown in FIG. 2, where represents 
the position of the quark qi (antiquark qC) which is de- 
noted by a solid (hollow) dot, the corresponding Jacobi 
coordinates can be expressed as 



r = r t - r 2 , R = r 3 - r 4 . 
miri + m 2 r 2 m 3 r 3 



X = 



TO4T4 



mi + m 2 



m 3 + mi 



(11) 



L, li and Z 2 are the orbital angular momenta associated 
with the relative motion coordinates X, r and R, respec- 
tively. The total wave function of a tetraquark state can 
be written as a sum of the following direct products of 
color, isospin, spin and spatial terms, 



[<^(r) Xsi ], tyg(R) X 



iJ 2 



■hi 



F l ( x )] j [VhVi 2 ]i [x Cl Xc 2 ] c , 



(12) 



where I and J are total isospin and angular momentum. 
a represents all possible intermediate quantum numbers 
a = {k,Si,Ji,Ji2,L,Ii}, where i = 1,2. x Si , m, an d 
Xa are spin, flavor and color wave functions of diquark 
or anti-diquark, respectively. [ ]'s denote Clebsh-Gordan 
coefficients coupling. The overall color singlet can be 
constructed in two ways, \\ = ^12 ® 3 3 <j, \\ = 6i 2 ® 6 3 4, 
"good" diquark and "bad" diquark are both included. 
Taking into account all degrees of freedom, the Pauli 
principle must be satisfied for each subsystem of identical 
quarks or antiquarks. The coefficient ^ is determined 
by diagonalizing the Hamiltonian. 

To obtain a reliable solution of a few-body problem, a 
high precision method is indispensable. In this work, a 
Gaussian Expansion Method (GEM) [6J|, which has been 
proven to be rather powerful in solving a few-body prob- 
lem, is used to study four-body systems in the flux-tube 
model. In GEM, three relative motion wave functions 
are expanded as, 



Y hmi (r), 
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C« 2 ( R ) = E c n2 N n2l2 R l ^e-^ R2 Y l2m2 {A){li) 



n 2 = l 



713 = 1 



where N ni i 1 , N ri2 i 2 and N n3 i 3 are normalization con- 
stants. Gaussian size parameters are taken as the fol- 
lowing geometric progression numbers 



r n =na n \ a 



(14) 



The parity for a diquark-antidiquark state is the prod- 
uct of the intrinsic parities of two quarks and two anti- 
quarks times the factors coming from the spherical har- 
monics 165] 



P = P q P q P q P q {-l) 1 ^ 



(-1) 



h+h+L 



(15) 



Using our coordinates, the eigenvalues of the charge con- 
jugation of a diquark-antidiquark state can be calculated 
by following the same steps as in the qq case. We can con- 
sider a diquark-antidiquark state as a QQ meson, where 
Q and Q represent a diquark and an antidiquark, respec- 
tively, with total "spin" J 12 and relative angular mo- 
mentum L between Q and Q (see Eq. 8). The C-parity 
eigenvectors are those states for which Q and Q have 
opposite charges. So applying the charge conjugation 
operator to these mesons is the same as exchanging the 
couple of quarks with the couple of antiquarks. The fac- 
tors arising from this exchange are the C-parity operator 
eigenvalues [§f|, 



C = (-l) 



L+J 12 



(16) 



The G-parity is a generalization of the concept of C- 
parity such that members of an isospin multiplet can 
each be assigned a good quantum number that would re- 
produce C-parity for the neutral particle. The G-parity 
operator is defined as the combination of C-parity and a 
7r rotation around the y axis in the isospin space [HJ , 



G=CKy{TT) =Ce 



(17) 



The G-parity eigenstates are tetraquark states with flavor 
charges equal to zero, i.e., strangeness equal to zero in the 
light mesons case, and their eigenvalues are: 



G=(-l) 



L+J 12 +I 



(18) 



V. NUMERICAL RESULTS AND DISCUSSIONS 

The diquark (antidiquark) is considered as a new com- 
pound object Q (Q) with no internal spatial excitations, 
and spatial excitations are assumed to occur only be- 
tween Q and Q in the present numerical calculations, 
which results in that such a tetraquark state has a lower 



energy than that of an internal spatial excited ones. The 
orbital angular momentum l\ and I2 are therefore as- 
sumed to be zero. With these restrictions the interme- 
diate quantum number J12 is the total spin angular mo- 
mentum S. The parity of a tetraquark with the diquark- 
antidiquark structure is P = (— 1) L , the charge conjuga- 
tion is C = (-1) L+S and the G parity is G = {-1) L + S+I . 

Within the flux-tube model with the parameters fixed 
by fitting the ordinary meson spectra [431 ] , the convergent 
energies of tetraquark states with this QCD quark cy- 
clobutadiene and diquark-antidiquark structures can be 
obtained by solving the four-body Schrodinger equation 



(H 



1 j 



0. 



(19) 



with Rayleigh-Ritz variational principle by setting the 
numbers of Gaussian wave functions n\ max = ri2max = 
n^max — 6. Minimum and maximum ranges of the bases 
are 0.1 fm and 2.0 fm for Jacobi coordinates r, R and 
X, respectively. Quark contents and the corresponding 
masses with specified quantum numbers, I G J PC or IJ P , 
are shown in Tables I-V, where n stands for a non-strange 
quark (u or d) while s stands for a strange quark, Ej and 
Ejj represent the energies of a QCD quark cyclobutadi- 
ene and a diquark-antidiquark structure, respectively, the 
quantum number ./V denotes the total radial excitation. 
The exotic in Table Til stands for a meson state which 
cannot be described by a qq configuration. 

The tetraquark states in the flux-tube model are gener- 
ally lower than that in the traditional quark models with 
additive two-body confinement interaction with color fac- 
tors used in early multiquark state calculations [43l |67| . 
The reason is that the multibody confinement potential 
can avoid the appearance of the anti-confinement in a 
color symmetric quark or antiquark pair. From Tables 
I-V, it can be seen that the two structures generally give 
very close energies for tetraquark ground states. How- 
ever, the differences between two structures are about 
40 MeV, 80 Mey and 120 MeV for spatial excitations 
between Q and Q with L = 1, L = 2 and L = 3, respec- 
tively, which is attributed to the reasons: (i) The expect 
values of the quadratic confinement potentials linearly 
depend on the angular excitation L [641 ] ; (ii) The normal 
modes of the three quadratic confinement potentials of 
two flux-tube structures, see Eq.(7) and Eq.(8) are dif- 
ferent. For a compact tetraquark state in ground state, 
the separation among particles (quarks or antiquarks) is 
generally smaller than one fm [42| , so the square of the 
length of each flux-tube is smaller than the length itself. 
It is therefore predicted that the small difference of the 
same quantum state between a linear confinement and a 
quadratic one is about 50-80 MeV according to the calcu- 
lations on hexaquark states 4l| . Anyway, the differences 
of the ground states between two structures are not big 
for both the linear and quadratic confinement potentials. 

In general, a tetraquark system should be the mix- 
ture of all possible flux-tube structures. Such as in the 
process of a meson-meson scattering, when two color sin- 
glet mesons are separated far away, the dominant com- 
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ponent of the system should be two isolated color singlet 
mesons because other hidden color flux-tube structures 
are suppressed due to the confinement. With the sep- 
aration reduction, a deuteron-like meson-meson molecu- 
lar state may be formed if the attractive force between 
two color singlet mesons is strong enough. When they 
are close enough to be within the range of confinement 
(about 1 fm), all possible flux-tube structures including 
the QCD quark cyclobutadiene and even more compli- 
cated flux-tube structures may appear due to the excita- 
tion and rearrangements of flux-tubes and junctions. All 
of these hidden color components can not directly decay 
into two colorful hadrons due to the color confinement. 
They must transform back into two color singlet mesons 
by means of the rupture and recombination of flux-tubes 
before decaying into two color singlet mesons. The de- 
cay widths of these states are qualitatively determined by 
the speed of the rupture and recombination of the flux- 
tubes. This formation and decay mechanisms are similar 
to the compound nucleus formation and therefore should 
induce a resonance which is named as a "color confined, 
multi-quark resonance" state (68| . It is different from 
all of those microscopic resonances discussed by S. Wein- 
berg [69( . Bicudo and Cardoso studied tetraquark states 
using the triple flip-flop potential including two meson- 
meson potentials and the tetraquark four-body potential. 
They also found plausible the existence of resonances in 
which the tetraquark component originated by a flip-flop 
potential is the dominant one (70j . 

Most tetraquark states in Tables I-V have the same 
quantum numbers with ordinary meson states, and the 
calculated energies of many tetraquark states are very 
close to the experimental data of the mesons with the 
same quantum numbers [7lj , especially states with higher 
energy. This does not mean that the main component of 
those experimental states must be tetraquark states. The 
fact is that most of the experimentally observed mesons 
can be interpreted as qq states (at least the main compo- 
nent) and accommodated in the naive quark model, only 
a few of them may go beyond qq configurations |72l l73j . 
However, the calculations indicate that the tetraquark 
component in those mesons, their energies are close to the 
tetraquark ones, can not be excluded. This point is sup- 
ported by the study on the nature of scalar mesons 74 1 . 
Moreover the nucleon spin structure study shows that 
even for the ground state the pentaquark components 
q qq is indispensable in solving the proton spin "cri- 
sis" [75], [76| . The strange magnetic momentum of a nu- 
cleon originates from a strange sea quark ss component 
is nonzero [77}. So a comprehensive study of the me- 
son spectra must include the mixing of qq and qqqq Fock 
components and in turn requires the knowledge of the 
off-shell interaction for annihilating or creating a quark- 
antiquark pair into or from the vacuum. Quark model 
should be unquenched and such an unquenched quark 
model study is on going in our group. 

With regard to nonstrange mesons, for some light qq 
excitation states, the orbital excitation energy between 



TABLE I: The mass spectra for nnnn states (unit:MeV). 



jGjPC 


N 2S+1 Lj 


Ej 


E ZI 


States 


PDG 


0+0++ 


l x 5 


601 


587 


/o(600) 


400-1200 


0+0++ 


2 1 S 


1101 


1019 


/o(980) 


980 ± 10 


0+1++ 


l 5 Di 


1927 


1840 


A(1285) 


1281.8 ± 0.6 


0+1++ 


2 5 Di 


1984 


1919 


A (1420) 


1426.4 ± 0.9 


0+1++ 


3 5 Di 


2373 


2270 


A (1510) 


1518 ±5 


0+2++ 


l 1 D 2 


1328 


1196 


A (1270) 


1275.1 ± 1.2 


0+2++ 


2 1 D 2 


1809 


1614 


A (1640) 


1639 ± 6 


0+2++ 


1 5 S2 


1468 


1465 


A (1430) 


w 1430 


0+2++ 


2 5 S 2 


1495 


1508 


A (1525) 


1525 ± 5 


0+2++ 


I s D 2 


1927 


1840 


A (1910) 


1903 ± 9 


0+2++ 


2 5 D 2 


1984 


1919 


A (1950) 


1944 ± 12 


O+4++ 


2 5 D 2 


1984 


1919 


A (2050) 


2018 ± 11 


(T2+- 


1 3 D 2 


1908 


1836 


exotic 


— 


0+0-+ 


l 3 Po 


1624 


1609 


77(1295) 


1294 ± 4 


o+o-+ 


2 3 P 


1656 


1619 


77(1405) 


1409.8 ± 2.5 


0+0-+ 


3 3 P 


2063 


2027 


77(1475) 


1476 ± 4 


0+0-+ 


4 3 P 


2097 


2055 


77(1760) 


1756 ± 9 


0+2-+ 


1 3 P 2 


1624 


1609 


772(1645) 


1617 ±5 


Q-l— 


l 1 ^ 


1057 


975 


0(1020) 


1019.455 ± 0.020 


0-1— 


2 1 P 1 


1482 


1358 


^(1420) 


1400-1450 


o-i~ 


3 x Pi 


1583 


1536 


w(1650) 


1670 ± 30 


0-1— 


l 5 Pi 


1696 


1651 


w(1650) 


1670 ± 30 


O-3— 


1 5 P 3 


1696 


1651 


u; 3 (1670) 


1667 ± 4 


o-i+- 


l 3 Si 


1291 


1304 


fci(1170) 


1170 ± 20 


0-1+- 


2 3 Si 


1391 


1394 


fci(1380) 


1386 ± 19 


1-0++ 


l^o 


1202 


1210 


a (980) 


980 ± 20 


1-0++ 


2 3 S 


1520 


1528 


a (1450) 


1474 ± 19 


1-1++ 


l 5 L>i 


1927 


1839 


ai(1260) 


1230 ± 40 


1-1++ 


2 5 Di 


2373 


2271 


ai(1640) 


1647 ± 22 


1-2++ 


1 5 S 2 


1470 


1467 


a 2 (1320) 


1318.3 ±0.6 


1-2++ 


l 1 D 2 


1876 


1807 


a 2 (1700) 


1732 ± 16 


1+2+- 


1 3 D 2 


1910 


1837 


exotic 


— 


i-o-+ 


l 3 Po 


1371 


1307 


tt(1300) 


1300 ± 100 


1-1-+ 


l 3 Pi 


1371 


1307 


7Tl(1400) 


1354 ± 25 


1-1-+ 


l 5 fi 


1775 


1691 


7Ti(1600) 


1662^f 


1+1— 


1 X A 


1580 


1558 


p(1570) 


1570 ± 36 ± 62 


1+1— 


1 5 F 1 


2157 


2030 


p(2150) 


2149 ± 17 


1+3— 


1 5 P 3 


1697 


1651 


p 3 (1690) 


1686 ± 4 


1+3— 


2 5 P 3 


2146 


2062 


p 3 (1990) 


1982 ± 14 


1+1+- 


l 3 Si 


1070 


1089 


6i(1235) 


1229.5 ± 3.2 


2+0++ 


l'So 


1202 


1211 


exotic 


— 


2+0- + 


l 3 Po 


1655 


1617 


exotic 


— 


2-1— 


l'Pi 


1580 


1558 


exotic 




2-l~ 


l 5 Pi 


1697 


1651 


exotic 




2-1+- 


l 3 Si 


1388 


1391 


exotic 




2+1++ 


l 5 Di 


1927 


1840 


exotic 




2+2++ 


1 1 D 2 


1876 


1807 


exotic 




2+2++ 


1 5 S 2 


1468 


1470 


exotic 
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q and q may be higher than that of a quark-antiquark 
pair excited from the quark sea, so these meson states 
prefer to have high Fock component qqqq. Such as the 
meson a, it can be described as the ground state with 
the quark content nnnn rather than the excited states 
of a qq meson 1431 , which is consistent with many other 
works [13, Fz8l - 480l j . The first radial excited state of the 
nnnn state is very close to the experimental value of the 
meson / (980), the tetraquark state nnnn may therefore 
be one of the main components, which is supported by 
Vijande's work on the nature of scalar mesons [7^ . The 
decay of the meson /o (980) into KK can be accounted for 
by other strangeness components, such as ss and nsns. 
The meson /o(1500) can not be described as a qq me- 
son, the mass and decay are compatible with it being 
the ground state glueball mixed with the nearby states 
of the ++ qq nonet [73|. In the quark models, another 
interpretation of the meson /o(1500) is that the main 
component might be a tetraquark state nsns [74j . The 
meson _£2(1430) has no proper member in the qq picture 
either [72 ]. It is suggested that the main component is 
a tetraquark nnnn, with quantum numbers 1 S2 in the 
flux-tube model. This state is not confirmed in the PDG 
and even recent measurements have suggested a differ- 
ent assignment of quantum numbers, which could make 
it compatible with the lightest scalar glueball [U . 

With respect to I = | strange mesons, most of them 
can be interpreted as the states which dominated by qq 
components in the quark models except for three mesons 
k(800), if*(1410) and if 2 (1580) [z3. For the same rea- 
son with the meson a, our model recommends a ground 
tetraquak state nnfis l l So with energy close to the meson 
k(800), which is compatible with other works [371. 1781480] . 
For the meson if*(1410), its assignment to the 2 3 5i state 
of the meson If* (892) is not only excluded by the large 
mass difference, but also by its decay modes [72J. A pos- 
sible interpretation of the main component of this state 
is a tetraquark state nnfis with quantum numbers 1 Pi 
instead of a pure qq pair. The meson ^(1580) has also 
no proper member in the qq spectra [72j . our tetraquark 
state l 3 .P2 mass is a little lower than experimental data. 
This state is clearly uncertain, it was reported in only 
one experimental work more than twenty years ago and 
has never been measured again. 

Concerning the exotic meson sector, the quantum num- 
bers rule out the pure qq possibility. The iri mesons of 

jGjPC = x - r 

~+are listed as manifestly exotic states 
by several experiments [82Tj85j . Many theoretical stud- 
ies have been made and various interpretations were 
proposed: hybrid meson states I86H89II. Trn molecular 
states [§(| and tetraquark states [9U |92|. Two mesons 
7Ti (1400) and tti (1600) are studied in the flux-tube model, 
see Table I, which indicates that the main components of 
7Ti(1400) and 7Ti(1600) might be tetraquark state nnnn 
with quantum numbers 1 3 P\ and 1 Pi, respectively. The 
tetraquark states nsns with quantum numbers 7 = 0,1 
and J PC = 1 h are predicted in the flux-tube model, 
the energies are around 1850 MeV (see Table II) which is 



TABLE II: The mass spectra for nsns states (unit:MeV). 



jGjPC 


N 2S+l Lj 


Ei 


En 


States 


PDG 


0+0++ 


^So 


1316 


1318 


/o(1370) 


1200-1500 


0+0++ 


2'So 


1583 


1590 


/o(1500) 


1505 ± 6 


0+0++ 


3 X S 


1676 


1661 


/o(mo) 


1720 ± 6 


0+0++ 


1 5 D 


2174 


2095 


/o(2100) 


2103 ± 8 


0+0++ 


l 5 D 


2174 


2095 


/o(2200) 


2189 ± 13 


0+2++ 


1 5 S 2 


1751 


1755 


/a(iaio) 


1815 ± 12 


0+2++ 




2033 


1946 


/a (2010) 


2011±?§ 


0+2++ 


2 1 Z> 2 


2141 


2073 


h (2150) 


2157 ± 12 


0+2++ 


1 5 D 2 


2174 


2095 


/a (2150) 


2157 ± 12 


o+o-+ 


1 3 P 


1867 


1831 


77(1760) 


1756 ± 9 


0+2-+ 


1 3 P 2 


1867 


1831 


t? 2 (1870) 


1842 ± 8 


0~1" 


l x Pi 


1773 


1740 


<^>(1680) 


1680 ± 20 


0~1" 


2 1 P 1 


1892 


1866 


— 


— 


o-i+- 




1583 


1586 


fci (1595) 


1594 ± 15^60 


o-i+- 


2 3 5'i 


1626 


1628 


— 


— 


o+i-+ 


l 3 Pi 


1865 


1828 


exotic 


— 


(T2+- 


1 3 Z> 2 


2108 


2076 


exotic 


— 


O-3— 


1 5 P 3 


1968 


1928 


03(1850) 


1854 ± 7 


1-0++ 


l'So 


1320 


1318 


a (980) 


980 ± 20 


1-0++ 


2'So 


1584 


1590 


a (1450) 


1474 ± 19 


1-2++ 


l 5 S 2 


1751 


1755 


a 2 (1700) 


1732 ± 16 


1-2++ 


1^2 


2033 


1945 


— 


— 


i-o-+ 


l 3 Po 


1867 


1831 


tt(1800) 


1816 ± 14 


1-1-+ 


l 3 Pi 


1867 


1831 


exotic 


— 


1-2+- 


1 3 D 2 


2108 


2076 


exotic 


— 


1-2-+ 


1 3 P 2 


1867 


1831 


tt 2 (1880) 


1895 ± 16 


1-2-+ 


1 3 P 2 


2309 


2186 


tt 2 (2100) 


2090 ± 29 


1+1 — 


l 1 ^ 


1772 


1739 


p(1700) 


1700 ± 20 


1+1 — 


2 1 P 1 


1892 


1866 


p(1900) 


1909 ± 17 ± 25 


1+1 — 


1 5 P 5 


2376 


2259 


p(2150) 


2149 ± 17 


1+3— 


1 5 P 3 


1967 


1928 


p 3 (1990) 


1982 ± 14 


1+3- 


l'Ps 


2248 


2117 


p 3 (2250) 


~ 2232 


1+5- 


1 5 P 5 


2376 


2259 


p 5 (2350) 


2330 ± 35 


1-1++ 


l 5 Di 


2173 


2095 






1+1+- 


l 3 ^! 


1583 


1586 







consistent with the predictions on J PC = 1 ^ tetraquark 
states in the QCD sum rule j9j, [12] ■ The exotic meson 
states with quantum numbers J PC = 2^ are predicted 
in the tetraquark picture, the states nnnn and nsns have 
the lowest masses around 1880 MeV and 2100 MeV, re- 
spectively. In addition, many flavor 1 — 2 exotic meson 
states are also calculated for further studies, see Table I. 

VI. SUMMARY 

The QCD quark cyclobutadiene, a new flux-tube struc- 
ture, is proposed in the framework of the flux-tube model. 
The flux-tube ring in the QCD quark cyclobutadiene can 
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TABLE III: The mass spectra for ssss states. 





at29-I-1 r 


£/ 


En 


States 


PDG 


0+0++ 




1919 


1925 


/o(2020) 


1992 ± 16 


0+0++ 


l 5 I>o 


2440 


2365 


jo (2330) 


2314 ± 25 


0+2++ 


1 5 5 2 


2051 


2044 


h (2010) 


20111^ 


0+2++ 




2423 


2354 


/ 2 (2300) 


2297 ± 28 


0+2++ 


1 5 D 2 


2440 


2365 


h (2300) 


2297 ± 28 


+ 2 ++ 


1 £>2 


2423 


2354 


/ 2 (2340) 


2340 ± 55 


0+2++ 


1 J I>2 


2440 


2365 


/ 2 (2340) 


2340 ± 55 


0+4++ 


1 5 Z>2 


2440 


2365 


/ 4 (2300) 


- 2314 


o _ i— 


l X Pl 


2201 


2176 


0(2170) 


2175 ± 15 


0+0-+ 


l 3 Po 


2232 


2195 


j?(2225) 


2226 ± 16 


o _ i— 


l 5 Pl 


2249 


2209 


0(2170) 


2175 ± 15 


o-i+- 


l 3 Di 


2432 


2359 







TABLE IV: The mass spectra for nnns states. 



IJ P 


N 2S+1 Lj 


Ei 


Eh 


States 


PDG 


2-0+ 


1^0 


995 


947 


#0(800) 


676 ± 40 


¥ + 


2 1 S 


1383 


1380 


#o*(1430) 


1425.6 ± 1.5 


¥ + 


l 5 A. 


2050 


1968 


#5(1950) 


1945 ± 10 ± 20 


±2 + 

2 Z 


1 5 D 2 


2050 


1968 


# 2 * (1980) 


1973 ± 8 ± 25 


¥ + 


1 5 D 4 


2050 


1968 


# 4 * (2045) 


2045 ± 9 




1 3 P> 


1514 


1451 


#(1460) 


~ 1460 


¥- 


2 3 P 


1739 


1697 


#(1630) 


1629 ± 7 




3 3 P 


1772 


1754 


#(1830) 


~ 1830 


2-1" 


l x Pi 


1430 


1367 


#*(1410) 


1414 ± 15 




2 1 Pi 


1709 


1666 


#*(1680) 


1717 ±27 


il + 
2 ± 


l 3 Si 


1254 


1233 


#i(1270) 


1272 ± 7 


il + 
2 x 


2 3 S'i 


1447 


1456 


#i(1400) 


1403 ± 7 


±1 + 
2 ± 


l 3 Di 


1749 


1644 


#i(1650) 


1650 ± 50 


±2+ 
2 Z 


1 5 S 2 


1603 


1601 


# 2 * (1430) 


1425.6 ± 1.5 


±2+ 

2 Z 




1685 


1573 


# 2 * (1430) 


1425.6 ± 1.5 


i2+ 

2 Z 


2 1 D 2 


2014 


1942 


# 2 * (1980) 


1973 ± 8 ± 25 


±2~ 
2 Z 


1 3 P 2 


1514 


1451 


# 2 (1580) 


- 1580 


±2" 
2 Z 


1 5 P 2 


1828 


1786 


#2 (1770) 


1773 ± 8 


i2" 

2 Z 


1 5 P 2 


1828 


1786 


#2 (1820) 


1816 ± 13 


2 J 


1 5 P 3 


1828 


1786 


#3 (1780) 


1776 ± 7 



be described as a glueball, four quarks are connected to 
the flux-tube ring by four fundamental flux-tubes, thus 
the QCD quark cyclobutadiene can be viewed as a qqqq- 
glueball hybrid. It provides a new sample for understand- 
ing the structures of exotic hadrons. The three familiar 
flux-tube structures can be taken as the ground states of 
a tetraquark system, the QCD quark cyclobutadiene may 
be an excited state which is obtained by means of creating 
Y-shaped junctions and flux-tubes from the vacuum and 
the rearrangement of some flux-tubes. The QCD quark 
cyclobutadiene and other three flux-tube structures are 
QCD isomeric compounds due to the same quark com- 
ponent and different flux-tube structures. 



TABLE V: The mass spectra for nsss states. 



ir 


N + L 


j Ei 


En 


States 


PDG 




l^o 


1757 


1762 


— 


— 




^ 1 s-< 

2 1 S 


1938 


1945 


#o(1950) 


1945 ± 10 ± 20 


1 ,-.-1- 


l°D 3 


2308 


2230 


# 3 (2320) 


2324 ± 24 




l 3 Po 


2026 


1984 


— 


— 


i „ 

2-0 


2 3 P 


2088 


2051 


— 


— 




l'Pi 


2051 


2024 








2 J Pi 


2160 


2139 








52- 


1 5 P 2 


2108 


2068 


# 2 * (2250) 


2247 ± 17 


2^ 


1 5 P 4 


2503 


2386 


#1(2500) 


2490 ± 20 


1 i — 


1 5 P 5 


2503 


2386 


#5(2380) 


2382 ± 14 ± 19 


±i+ 

2 ± 


I 3 Si 


1778 


1774 






±1+ 
2 ± 


2 3 S' 1 


1864 


1862 






\2+ 

2 




1904 


1900 






±2+ 

2 Z 


l 1 D 2 


2284 


2215 






±2+ 


2 l D 2 


2440 


2386 






Most 


meson 


states 


in the PDG can 


be described 



qq configurations, only a few meson states, a, k(800), 
/o(980), / (1500), 7ri(1400), 7ri(1600), / 2 (1430) and 
#*(1410), are difficult to be interpreted as qq mesons. 
The tetraquark state as their main components is one of 
possible interpretations of their flavor components. Some 
exotic meson states as tetraquark states are predicted in 
the flux-tube model. The tetraquark states nsfis with 
quantum numbers I = 0, 1 and J PC = 1 h have the 
lowest masses around 1850 MeV. The tetraquark states 
nnnfi and nsns with quantum numbers J PC = 2 H have 
the lowest masses around 1880 MeV and 2100 MeV, re- 
spectively. The tetraquark states with 1 = 2 are also 
predicted in the flux-tube model. 

Even though up to now no tetraquark state has been 
well established experimentally. It is indispensable to 
continue the study of the tetraquark system because the 
tetraquark component in mesons can not be ruled out 
and may play an important role in the properties of 
mesons, similar to the fact that the pentaquark compo- 
nents play an important role even in the nucleon ground 
state. 

The tetraquark states, if they really exist, should be 
the mixtures of all kinds of flux-tube structures which can 
transform one another. In this way, the flip-flop of flux- 
tube structures can induce a resonance which is named as 
a "color confined, multi-quark resonance" state. To ver- 
ify such a new resonance is not easy. We admit that this 
analysis is based on the mass calculation only, the crucial 
test of the components of exotic mesons is determined by 
the systematic study of their decays, which involves chan- 
nel coupling calculation containing all possible flux-tube 
structures and mixing between qq and tetraquark compo- 
nents and so much more information of low energy QCD 
are needed. 
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